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Abstract Beside its intrinsic interest for the insights it can give into color confinement, knowledge of the 
space-time evolution of hadronization is very important for correctly interpreting jet-quenching data in 
heavy ion collisions and extracting the properties of the produced medium. On the experimental side, the 
cleanest environment to study the space-time evolution of hadronization is semi-inclusive Deeply Inelastic 
Scattering on nuclear targets. On the theoretical side, 2 frameworks are presently competing to explain the 
observed attenuation of hadron production: quark energy loss (with hadron formation outside the nucleus) 
and nuclear absorption (with hadronization starting inside the nucleus). I discuss recent observables and 
ideas which will help to distinguish these 2 mechanisms and to measure the time scales of the hadronization 
process. 
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1 Introduction 

One of the most striking experimental discoveries in the 
heavy-ion program at the Relativistic Heavy Ion Collider 
(RHIC) has been the the suppression of large transverse 
momentum hadron production in nucleus-nucleus (A-l-A) 
collisions compared to proton-proton (p+p) collisions 1 . 
The observable of interest is the ratio of the hadron trans- 
verse momentum (pt) spectrum in A-l-B collision in a 
given centrality class (cel.), normalized to binary scaled 
p-|-p collisions by the inverse thickness function Tab and 
finally divided by the p-|-p spectrum: 



Rab = 



1 



Tab{c.cI. 



^^A+B^h+X 
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da 



p+p^h+X 



dp^dy 



dp:^dy 
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In the absence of nuclear effects, one would expect Rab — 
1. By comparing the measured RauAu ~ 0.2 in Au-|-Au 
collisions to the mild deformation of pr spectra observed 
in deuteron-gold (d-|-Au) collisions, one concludes that the 
large suppression of RauAu is due to the hot and dense 
medium produced in Au+Au collision, also called "hot 
nuclear matter" , see Fig.^ This measurement is one of the 
keys to the claimed discovery of the Quark-Gluon Plasma 
(QGP) at RHIC 2 . 

The suppression of hadron production in A+A colli- 
sions, has been successfully described in terms of parton 
energy loss due to medium-induced gluon bremsstrahlung, 
allowing so-called "jet tomography" studies of the QGP 
However, this success has been recently questioned. 
Gluon radiation off heavy quarks is expected on theoret- 
ical grounds to be suppressed at small angles compared 
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to light quarks, implying a smaller suppression for D and 
B mesons than for tt mesons However, the measured 
suppression of single non-photonic electrons at RHIC 
which are the decay product of B and D mesons, is of 
similar magnitude for pions contrary to theoretical expec- 
tations [7] . The common assumption of neglecting elastic 
parton energy loss in Raa computations has been recently 
reexamined j8, , but the effect seems insufficient to explain 
the data, at least within conventional schemes for treating 
the running coupling constant j^. As further assumption 
that needs to be tested and will be addressed in this pa- 
per, is that the quark which fragments into the observed 
hadron traverses the whole medium and hadronizes well 
outside it. If untrue, in-medium interactions and screening 
of the hadronizing system would need to be accounted for. 

Hadron suppression has also been observed in fixed 
target Deeply Inelastic lepton- nucleus Scattering (nDIS). 
In this case, the medium which induces the attenuation 
is the target nucleus itself, also called "cold nuclear mat- 
ter", see FigQ] Experimental data are usually presented 
in terms of the "multiplicity ratio" [Tr)IITl.. l2 UA^ 



RU^h) 



1 d N'Xjzh) 
dzh 



dN^izh) 



dzh 



(2) 



i.e., the single hadron multiplicity on a target of mass 
number A normalized to the multiplicity on a deuteron 
target as a function of the hadron's fractional energy Zh = 
Efi/v, where v is the virtual photon energy. The ratios in 
the numerator and denominator cancel to a large extent 
initial state effects like the modifications of parton dis- 
tribution functions due to shadowing and EMC effects, 
exposing the nuclear modifications of the fragmentation 
process: if we assume factorization formulae to be valid. 
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Figure 1. Left: Lepton- nucleus scattering. The hadronizing 
quark travels through the target nucleus. Right: Nucleus- 
nucleus scattering. The hadronizing parton travels through the 
medium produced in the collisions. 



we have at leading order Rm ~ ^h/Dj^, i.e., the ratio 
of the fragmentation functions (FF) in the nucleus A and 
in deuteron. If no nuclear effects modify the fragmenta- 
tion process we would expect Rm ~ 1. In fact, what 
is experimentally observed is hadron suppression in the 
Zh = 0.2 - 1 and ly = 2 ~ 20 GeV range at HERMES 
[TTltnitn] . and in the = 20 - 200 GeV at the EMC ex- 
periment 201 • The flavor dependence of the multiplicity 
ratio has also been measured [T^ltOj , showing suppression 
for pions, kaons and antiprotons. Protons are enhanced at 
Zh ^ 0.4 and suppressed above: this is a "proton anomaly" 
analogous to the "baryon anomaly" observed in p-f-A and 
A-l-A collisions [T1ll5| . Both the quenching and the en- 
hancement increase with A. Data binned in v, and in the 
photon virtuality are also available from HERMES. 
Very high-statistics measurements will be available in the 
near future from the CLAS experiment at Jefferson Labs 
dH, with some preliminary results already presented [TB^i . 

The cleanest environment to address nuclear modifi- 
cations of hadron production is nuclear DIS: it allows to 
experimentally control many kinematic variables; the nu- 
clear medium (i.e., the nucleus itself) is well known; the 
multiplicity in the final state is low, allowing for precise 
measurements. Moreover, the nucleons act as femtometer- 
scale detectors of the hadronizing quark, allowing to ex- 
perimentally study its space-time evolution into the ob- 
served hadron. Hadron suppression at HERMES and CLAS 
is of direct relevance to RHIC physics. In both cases the 
hadronizing quark has to traverse a length of (hot or cold) 
nuclear matter of the size of a nuclear radius, see Fig^ 
Moreover, in the HERMES experiment « 2 - 25 GeV 
and z^i w 0.2 — 1, so that hadrons are observed over an 
energy range = Zhi' = 2 — 20 GeV. (measurements at 
CLAS with Ebeam = 5, 11 GeV will likewise explore the 
Eh — 2 — IQ GeV range). For midrapidity hadron produc- 
tion in A+A coUisions at RHIC Eh k. pT k, 2 - 2Q GeV, 
roughly equal to the hadron energy range at HERMES. 

Information about parton propagation in cold nuclear 
matter is needed as an input for the interpretation of data 
in A-l-A collisions. In this case one wants to use hadron 
suppression as a tool to extract the properties of the hot 
QGP created in the collision. To this purpose we need to 
develop well calibrated computational tools to relate the 



magnitude of hadron suppression to properties of the QGP 
like its density and temperature. Assuming long lived par- 
tons, hadron suppression at RHIC would be attributed 
to parton energy loss, leading to a medium temperature 
of T w 400 MeV jlj, in excess of the critical tempera- 
ture Tc « 170 MeV for the QGP phase transition. If, on 
the contrary, hadronization started on the nuclear radius 
scale or before, in-medium interactions should also be ac- 
counted for, leading to a different, presumably lower T . A 
precise knowledge of parton propagation and hadroniza- 
tion mechanisms obtained from nDIS data is essential for 
testing and calibrating our theoretical tools, and to deter- 
mine the properties of the QGP produced at RHIC. 



2 Formation time estimates 

The key quantity we need to investigate is the hadroniza- 
tion time scale. Since hadronization is a non perturbative 
process, one has to resort to phenomenological models to 
describe it |17lll8lll9ll20ll21ll22ll2:^j . However, a few features 
are expected on general grounds. Due to color confine- 
ment, the struck quark must neutralize its color at some 
stage, say by picking up an antiquark from the vacuum or 
the surrounding medium. I call this color neutral qq pair a 
"prehadron" /i,, and the time for its formation the "pre- 
hadron formation time" (some authors prefer to call it 
the "production" time). This is a relevant time scale since 
gluon bremsstrahlung off the struck quark stops after color 
neutralization; moreover, the prehadron quickly develops 
a cross section of the order of the hadronic one, leading to 
its nuclear absorption. Subsequently, the prehadron wave 
function collapses on the observed hadron h wave function, 
and the corresponding time is called "hadron formation 
time" th- Pi. final caveat: it is difficult to rigorously define 
the concept of formation time in field theory, so that in 
the following discussions it is used as a working tool. 



2.1 Long formation time: energy loss models 

The average hadron formation time {th) can be considered 
as the time for the struck partons to build up its color 
field and to develop the hadronic wave function [^J. In 
the hadron rest frame this time is related to the hadron 
radius, and in the laboratory frame it is boosted to: 



Eh zv 

{th) oc Rh = Rh 

TUh mh 



(3) 



For a 10 GeV pion at HERMES, we have (th) ~ 50 fm 
^ Ra- Note also that the scale for hadron formation is set 
by Kh — mh/Rh ~ 0.2 GeV/fm. This estimate is used in 
energy loss models [TTItT^ to assume that hadronization 
starts outside the medium with a decreased parton en- 
ergy due to multiple parton scatterings and induced gluon 
bremsstrahlung. These models are fairly successful in de- 
scribing Rm at HERMES, see Fig.|21 

More in detail, the Ref. ^7] computes parton rescat- 
terings and gluon radiation in pQCD including Feynman 
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Figure 2. Left: Multiplicity ratio for at HERMES |12) compared to the energy loss model of Ref. |17| and the absorption 
models of Ref. |21I23| . Center: pure absorption model |21I27| (solid) and energy loss model .27I18| (dashed). Right: Color dipole 
model from Ref. |22| (dashed: absorption only, solid: absorption and induced energy loss). Data are from Ref. |11I12I13| . 



diagrams up to twist-4. Fragmentation of both the struck 
quark and the radiated gluon is included. The modifica- 
tion of the FF depends on 1 parameter, the strength of 
parton-parton correlations in the nucleus. The modified 
FF so obtained can be modeled to a good accuracy by 
shifting z/i in the leading-twist fragmentation function: 



1 



-D 



l-Azh V 1 - Azh 



Zh 



(4) 



where Azh = AEq/v is the quark's fractional energy loss, 
and AEq m 0.6(zg) with {zg) the average fractional energy 
carried away by the radiated gluon. 

In Ref. TS" the parton energy loss is treated in the 
BDMS formalism 25 , which takes into account medium- 
induced multiple soft gluon emission and fluctuations in 
the energy loss. Modified FF are computed as the average 
of the Zft-shifted FF in Eq. Q: 



(l-z) 



dAz P{Az;q, Lg 



I- A: 



-D){ 



Zh 



1- Az 



The "quenching weight" V{Az) is the probability distri- 
bution of a fractional energy loss Az = AEq/v |2.5I.!26? , 
and Lq the quark's in-medium path length. The medium 
is characterized by the transport coefficient g, which mea- 
sures the average momentum transfer per unit path-length 
from the medium to the parton. When also taking into ac- 
count a realistic medium geometry and finite medium size 
corrections to the quenching weight, the model ,27,can 
well describe HERMES data, see Fig. 12 



2.2 Short formation time: nuclear absorption models 

A successful non perturbative model of hadronization is 
the Lund string model 28 . The confined color field stretch- 
ing from the struck quark to the rest of the nucleus is 
modeled as a string of tension Kstr ~ 1 GcV/fm. Pre- 
hadron formation is identified with the qq pair production 



which breaks the string in smaller pieces ■ Hadrons are 
formed when a quark and an antiquark at the endpoint of 
a string fragment meet. Average formation times can be 
analytically computed [roil2UU21| : 



ZhV 
^str 



{th) = (i*) + 



^str 



(5) 
(6) 



The factor Zhi' can be understood as a Lorentz boost fac- 
tor. The (1 — Zh) factor is due to energy conservation: 
a high-z;i hadron carries away an energy Zhi'; the string 
remainder has a small energy e = (1 — Zh)v and cannot 
stretch farther than L — e/nstr- Thus the string break- 
ing must occur on a time scale proportional to 1 — Zh- The 
function f{zh) is only a small deformation of (t,). At HER- 
MES {U) « 4 fm < i?A and {th) 6 - 10 fm > Ra- The 
hadron is typically formed at the periphery or outside the 
nucleus so that its interaction with the medium is negligi- 
ble. However, the prehadron is formed well inside and can 
start interacting with the nucleus. The nuclear absorption 
model of Ref. 21 can successfully explain Rm measure- 
ments at HERMES in terms of prehadron-nucleon inelas- 
tic scatterings with the above formation times estimate, 
see Fig. |5| The prehadron-nucleon inelastic cross-section 
is (T,(i^) — 0.35(T/i(i^) proportional to the experimental 
hadron- nucleon cross section ah ■ The proportionality fac- 
tor is fitted to 7r+ production data on a Kr target ■ The 
prehadron survival probability S*, is computed in terms of 
transport equations. Neglecting hadron absorption. 



d^bdy pA{b,y) x dx 



f dsApA {b,s) 



where (6, y) is the 7*-q interaction point, pA is the nuclear 
density, and one can recognize exponential probability dis- 
tributions for prehadron and hadron formation. 

In Ref. 122] the formation of a leading hadron {zh > 
0.5) is described in a pQCD inspired model. The struck 
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Figure 3. Left: results of the Rm = cA" fit for {E[e,N,Ne,Kr} 
at z=0.65 (solid: absorption; dashed: energy loss; dotted: 
data |12p . Right: computations including only heavy nuclei 
{Kr,Sn,Xe,W,Au,Pb}. 

quark radiates a gluon according to the Bethe-Heitler ra- 
diation spectrum. The gluon then splits into a qq pair, 
and the q recombines with the struck q to form the lead- 
ing prehadron. Medium interaction and evolution of the 
prehadron into the observed hadron is computed in a path- 
integral formalism for dipole propagation. The prehadron 
formation time is identified with the time at which the 
gluon becomes decoherent with the struck quark. The 
probability distribution in the prehadron formation time 
can be computed, and the average (t*) is 

{U)^{l-z,)^. (7) 

The scale is set by Kdip = « 10 GeV/fm at HERMES, 
and (i,) < 5 fm at Zh > 0.5. This model can successfully 
describe leading hadron suppression, see Fig. [fright. 

In Ref. |23| a different space-time picture of hadroniza- 
tion is advocated. Prehadrons are formed at t.^, = 0, and 
hadrons are formed at = {Eh/mh)TQ with tq = 0.5 fm. 
The leading prehadron interacts with the medium with 
a reduced hadronic cross-section. Subleading prehadrons 
do not interact with the medium until hadron formation. 
This picture is then embedded in a Monte Carlo trans- 
port model. A good description of HERMES data can be 
achieved, see Fig. |21 

3 Can we distinguish energy loss from hadron 
absorption? 

Most of the difference in the time estimates quoted in the 
previous Section lies in the different scale k used. E.g., 
Kh ~ 0.2Kst,. leads to the rather large (th) « 50 fm quoted 
in the energy loss model estimate instead of (th) ~ 10 fm 
quoted in the Lund model estimate. In the second case 
there would be no justification for neglecting the interac- 
tions of the forming hadron field with the nucleus. As the 
choice of the scale k is a debatable and model-dependent 
matter |22II29| . it is very important to look for observ- 
ables which are able to distinguish energy loss models 
and absorption models, or to directly detect in-medium 
hadronization effects. 



3.1 Mass number dependence 

In first approximation, one expects 1 — R'^j oc A^/"^ in en- 
ergy loss models because the average energy loss AEq oc 
(Lq) oc A^/^, due to the Landau-Pomeranchuk-Migdal in- 
terference in QCD On the other hand, in absorp- 
tion models the survival probability is proportional to the 
amount of traversed matter, so that 1 — i?^^^ oc {La) oc 
A^^^. Therefore a simple analysis of the A-dependence of 
should clearly signal which model is correct. 

This argument fails for absorption models |21II30| . If 
the prehadron were produced always at the 7*-quark inter- 
action point (i.e., = 0) then Rm = cA^^^ at all orders in 
A^^^ . However, if we allow for a nonzero (i*), its dimension 
must be neutralized by the nuclear radius Ra, introduc- 
ing extra powers of A^'^. Quite generally, if the probability 
distribution for the prehadron formation length is finite at 
zero formation length, then oc A'^^^ + 0{A), the same 
power found in energy loss models. 

Then, we can study the breaking of the A'^^^ law. To 
this purpose, it was proposed in ^TT to select a set of tar- 
gets {Ai,A2, ■ ■ ■ ,An}, fix the z bin, and perform a fit of 
the form 1 — R\j{z) — c{z)A°'^^'> with c and a free pa- 
rameters. Results are presented in terms of 2a confidence 
contours in the (c, a) plane. As shown in Fig.|31left, energy 
loss 127] and absorption models |21II27| are indistinguish- 
able. The same holds true for all z bins. Experimental data 
are described by an A"' power law with a = 0.61 ± 0.14, 
compatible with a = 2/3 but excluding a = 1/3. Increas- 
ing the number of targets and the span in atomic num- 
ber does not help in separating the 2 models, either, but 
clearly shows a non negligible breaking of the A^^"^ law at 
A > 80 21,27,, see Fig.l^right. 

3.2 Formation time scaling 

In Ref. 1^ I conjecture that Rm should not depend on 
Zh and V separately but should depend on a combination 
of them: 

Rm = RM[T{zh,iy)] , (8) 

where the scaling variable r is defined as 

T = Cz^(l-z,Oz^ . (9) 

The scaling exponent A is introduced as a way of approx- 
imating and summarizing the scaling behavior of exper- 
imental data and theoretical models. It can be obtained 
by a best fit analysis of data or theoretical computations. 
The proportionality constant C cannot be determined by 
the fit. A possible scaling of Rm with is not consid- 
ered in this analysis because of its model dependence, and 
because of the mild dependence of HERMES data on Q^. 
As discussed below, the proposed functional form of r, 
Eq. (O, is flexible enough to encompass both absorption 
models and energy loss models. The 2 classes of models 
are distinguished by the value of the scaling exponent: a 
positive A ^ is characteristic of absorption models, while 
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Figure 4. Left: An example of the fit procedure described in Section 13.21 applied to HERMES data for vr"*" production on a 
Kr target |12| . Upper panel: as a function of A. Lower panel: Rm{t) with r computed at Abcst ~ 0.34. Right: The scaling 
exponent Abest extracted from HERMES data on charged and identified hadrons at Eiab ~ 27 GeV |ilk.l2 13 (only statistical 
errors included in the fit). Error bars correspond to 1 standard deviation. The bottom panel shows the per degree of freedom. 



a negative A < is characteristic of energy loss models. 
Thus, the exponent A extracted from experimental data 
can identify in the experimental data the leading mecha- 
nism for hadron suppression in nDIS. 

The scaling of Rm is quite natural in the context of 
hadron absorption models |19II2()II2 1 22 23 . Indeed, pre- 
hadron absorption depends on the in-medium prehadron 
path length, which depends solely on the prehadron for- 
mation time (i*) as long as (i*) < Ra- As argued in Sec- 
tion |21 (i*) oc f{zh){l — ZhjZhV which is well described by 
the proposed scaling variable t with A > 0. E.g., in the 
Lund model A « 0.7. 

In energy loss models |17lll8ll2'7j the scaling is less obvi- 
ous. For the purpose of discussing the scaling properties of 
Rm , we can consider the model of Ref . ^18„.27j and neglect 
finite medium size corrections and finite quark energy cor- 
rections. If we further neglect energy loss fluctuations, we 
can approximate Rm ~ DA{zh)/ D(zh) and obtain 



Rm 



where the average energy loss (e) = /^^ ^''^"^ de e 'P(e) 

/ /o^ ^''^'^ deV{e) — f[{l~Zh)iy] is a function of the energy 
(1 — Zh)v not carried away by the observed hadron. Next, 
we can approximate the FF using the parametrization of 
Ref. at Q2 = 2 GeV^ and obtain 



Rm 



1 



^ 1 \ a+/5+l 

1 - -/[(I - z,);.] 
^ 1/ 



1 - 



/[(I - ^/OH 

(1 - Zh)!^ 



This shows an approximate scaling with (1 — z^)!^, which 
implies scaling of Rm with respect to r with A w 0. A 
similar argument holds for the model of |17| . When per- 
forming the scaling analysis of the full models one finds in 
general A < pTT] . 



The HERMES experiment measures Rm binned in 
Zh and integrated over v and ("z/j distributions") or 
binned in v and integrated over Zh and ('V distri- 
butions"). Eq. is fitted to the combined Zh- and i'- 
distributions, and the scaling exponent A is determined 
by minimization. An example of this procedure is il- 
lustrated in Fig. ??. Details can be found in Ref. (31) . 

The scaling exponents Abest extracted from HERMES 
data at Eiab = 27 GeV [TTlllijlll^l for different hadron 
flavors produced on N, Ne and Kr targets are shown in 
Fig. ??. In all cases x^/d.o.f. < 1.6, which proves that 
Rm scales with r. The central result of this analysis is 
that pion data exhibit a clear Abest ~ 0.4 ^ 0. As dis- 
cussed, this shows the dominance of the prehadron ab- 
sorption mechanism as opposed to the energy loss mecha- 
nism, or in other words is a signal of in- medium prehadron 
formation, with formation times (t^) < Ra- 



3.3 pt broadening 

The scaling analysis just described gives only indirect ev- 
idence for a short formation time, and cannot measure 
its absolute scale. An observable which is more directly 
related to the prehadron formation time is the hadron's 
transverse momentum broadening in DIS on a nuclear tar- 
get compared to a proton or deuteron target 22- Indeed, 
when a hadron is observed in the final state, neither the 
quark nor the prehadron could have had inelastic scatter- 
ings. The prehadron- nucleon elastic cross section is very 
small compared to the quark cross section, so that the 
hadron's p^-broadening originates doniinantly during par- 
ton propagation. As shown in (33^34 . the quark's momen- 
tum broadening Ap"^ is proportional to the quark path- 
length in the nucleus. If the prehadron formation time has 
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the form © as argued in the last section, we obtain: 
Ap\ oc (t,) oc z,-^(l - Zh)v , 

unless the distance of the quark production point from the 
surface of the nucleus is smaller than (i*). Then we should 
expect a decrease of Ap^ with increasing z^- This would 
be a clear and model-independent signal of in-medium 
prehadron formation: indeed, if the quark were travel- 
ing through the whole nucleus before prehadron forma- 
tion Ap^ would only depend on the nucleus size and not 
on Zh- A related observable is the z^-dependence of the 
Cronin effect, which is likewise expected to decrease with 
increasing Zh |22j . 

The CLAS collaboration can perform multi-differential 
PT-broadening measurements in all kinematic variables 
thanks to a very high beam luminosity. A few prelimi- 
nary data from CLAS are already available U^I ■ The HER- 
MES collaboration is also studying the pT-broadening at 
a larger beam energy but with a lower statistics. The scal- 
ing analysis proposed in the previous section will be useful 
to cross-check the results of these measurements. 



4 Conclusions and perspectives 

Use of hard processes to probe medium processes in A-l-A 
collisions requires a detailed understanding of the hadroniza- 
tion process, which can be studied in lepton-nucleus scat- 
terings. The scaling analysis 31 of pion attenuation at 
HERMES demonstrates for the first time a scaling of the 
hadron attenuation ratio Rm which is compatible with a 
short prehadron formation time of the order or smaller 
than the nuclear radius. Thus, it favors prehadron ab- 
sorption as dominant mechanism for hadron suppression 
instead of gluon radiation off a struck quark. This conclu- 
sion will be soon checked by hadron py-broadening data. 

Much more can be studied in lepton-nucleus scatter- 
ings, (i) In the meson sector, the suppression of 77 mesons 
at RHIC, which is of similar magnitude than for tt, seems 
to favor long-lived partons 35j. Measuring ry and heav- 
ier meson attenuation at HERMES and CLAS will check 
the correctness of such interpretation, (ii) Understanding 
the proton anomaly in nDIS will shed light on baryon 
transport in nuclear matter and on the baryon anomaly 
observed in p-f A and A-t-A collisions. Measurements of A 
and other baryons at HERMES and CLAS will be needed 
in this respect, (iii) Neither HERMES nor CLAS are able 
to study the hadronization of heavy quarks, because of 
limited luminosity and limited Bjorken's x coverage, re- 
spectively. The proposed Electron-Ion Collider 36 is well 
suited for such studies, thanks to its low- a; coverage and 
high luminosity. Study of D and B meson suppression will 
settle the single electron puzzle at RHIC and will put in- 
terpretation of LHC data on a firmer ground. Study of 
"normal" J /tp suppression will help in distinguishing com- 
peting mechanism and in building a precise baseline for 
measurements of the "anomalous" suppression in A+A 
collisions. 
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